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Abstract. We determine conditions under which a generic gauge invariant nonautonomous and 
inhomogeneous nonlinear partial differential equation in the two-dimensional space-time continuum 
^ ■ can be transform into standard autonomous forms. In addition to the nonlinear Schrodinger equa- 

tion, important examples include the derivative nonlinear Schrodinger equation, the quintic complex 
■ Ginzburg-Landau equation, and the Gerdjikov-Ivanov equation. This approach provides a mathe- 

matical description of nonstationary media supporting unidimensional signal propagation and/or to- 
tal field trapping. In particular, we study self-similar and nonspreading wave packets for Schrodinger 
\ equations. Some other important coherent structures are also analyzed and applications to nonlin- 

ear waves in inhomogeneous media such as atmospheric plasma, fiber optics, hydrodynamics, and 
fH ! Bose-Einstein condensation are discussed. 
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1. Introduction 



The nonlinear Schrodinger equation is a universal mathematical model for many physical systems 
— in a variety of nonlinear problems, the one- dimensional nonlinear Schrodinger equation with 
^ ', a cubic nonlinearity describes the propagation of an envelope wave riding over a carrier when 
• ; the lowest order contribution to the dispersion relation is proportional to the square of a local 
^ ■ amplitude of the wave envelope. This equation is generic arising in consideration of the lowest 
CN . effects of dispersion and nonlinearity on a wave packet, including nonlinear optics, deep water 
^ [ waves, acoustics, plasma physics and nonlinear propagation of plane diffracted wave beams in the 
focusing regions of the ionosphere [1], [IB], [IB], [27], [25], [71], [M], [HS], [S], [SH], [ISS], [M], [HZ], 
><■ [203]. and [204j. 

d \ The inverse scattering method is a standard approach to several completely integrable nonlinear 
partial differential equations, including the Korteweg-de Vries, nonlinear Schrodinger or Gross- 
Pitaevskii, Sine-Gordon and Kadomtsev-Petviashvili equations [2], [S], |143] . |154] . |165] . and |166] . 
At the same time, the physical situations in which these equations arise tend to be highly idealized. 
The inclusion of effects such as damping, external forces, an inhomogeneous medium with variable 
density and a higher order of the nonlinearity may provide a more realistic model. However, the 
addition of these perturbation effects could mean that the system is no longer completely integrable 
(see an example in Refs. [2], [31], [61], and |180j ). Hence it is of interest to determine under what 
conditions the perturbation preserves the completely integrability. Some nonintegrable systems 
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in question posses important in practice classes of exact solutions (see, for instance, |39], |I0] . 
[128] . |137] and the references therein). These solutions may serve as a natural starting point of 
perturbation methods and also provide a useful testing ground for numerical investigation of more 
complicated models. 

In this paper, we show that the following nonlinear PDE in the two-dimensional space-time 
continuum: 



where Q is a quadratic of two (noncommuting) operators d/dx and x with time-dependent coeffi- 
cients, the asterisk denotes the complex conjugation, and P is certain gauge invariant nonlinearity, 
can be reduced by a simple change of variables to the autonomous form: 



which is somehow easier to analyze by a variety of available tools (in general, the new time variable 
r may be complex-valued and the new space variable ^ is a linear complex-valued function of x; the 
precise form of this transformation is given by Theorem 1 below). Among important in applications 
examples are the derivative nonlinear Schrodinger equation, which describes the propagation of 
circular polarized nonlinear Alfven waves in plasma physics, the quintic complex Ginzburg-Landau 
equation, the Gerdjikov-Ivanov equation, and others. Here, we mainly concentrate on variants of 
the nonlinear Schrodinger equation for which the unperturbed models are known to be completely 
integrable and/or have explicit solutions. 

The paper is organized as follows. In the next section, we transform (11. ip into the standard form 
(II. 2p . An overview of integrable autonomous cases and some exact solutions of important nonin- 
tegrable PDE models are given in sections 3 to 6 in order to make our presentation as self-contain 
as possible (a detailed bibliography is provided). Some applications are discussed in section 7. In 
particular, we study self-similar and nonspreading wave packets for cubic nonlinear Schrodinger 
equations and quintic complex Ginzburg-Landau equations when some other interesting coherent 
structures are also available. A self-focusing wave packet for the linear Schrodinger equation, an 
"Airy gun", which self-accelerates to an infinite velocity in a finite time, is found in analogy to 
the well-known blow up property of nonlinear PDE solutions. A natural hierarchy of these wave 
packets is discussed from a "hidden" symmetry view point and applications to inhomogeneous me- 
dia including fiber optics, atmospheric plasma, hydrodynamics, and Bose condensation are briefiy 
reviewed. 

2. Transforming Nonlinear Schrodinger Equation into Autonomous Forms 
We start from the nonautonomous nonlinear Schrodinger equation 



on M, where the variable Hamiltonian H = Q {p, x, t) is an arbitrary quadratic form of two operators 
p = —id/dx and x, namely. 



= -a (t) + b (t) - it) x^, -id{t)ij-f [t) xi) + %g (t) + h (t) \i;f i, (2.2) 








(2.1) 
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(a, b, c, d, f, and g are suitable real- valued functions of time only) under the following integrability 
condition [T7i]. fm\: 

h = hoa (t) it) (t) = hoP' (0) ^2 (0) -^^J^^l^ (2.3) 

{ho is a constant, functions /3, A, and /i will be defined below). 

(The corresponding linear quantum systems, when h = 0, are known as the generalized (driven) 
harmonic oscillators. Some examples, a general approach and known elementary solutions can be 
found in Refs. |33], IS], [51], [123], USE], [M], and [2115] .) 

For completeness we combine the results established in |116] and |180] . 
Lemma 1. T/ie substitution 

i[a{t)x^+5(t)x+K(t)) 

^ = ^== X(e,r), e = /3(t)a: + £(t), r = 7(t) (2.4) 

> (t) 



transforms the non- autonomous and inhomogeneous Schrddinger equation into the autonomous 
form 

iXr + ho IxT X = Xi^- co^^X (co = 0, 1) (2.5) 

provided that 

— + b + 2ca + Aaa^ = (2.6) 

^ + (c + 4aa)/3 = 0, (2.7) 



dt 

+ a/32 = (^2.8) 



and 



Here 



dt 
dS 

— + (c + 4aa) S = f + 2ga + 2coa(3^e, (2.9) 

= (g-2a6)P, (2.10) 

at 

— = gS -a6'^ + coa/3^e'^. (2.11) 

1 a' d , , 

a = — — . 2.12 

4a 2a 



(A Mathematica based proof of Lemma 1 is given by Christoph Koutschan |108j .) 

The substitution (12.121) reduces the inhomogeneous equation (12. 6p to the second order ordinary 
differential equation: 

fx" - T (t) fx' + Aa (t) fx = Co {2af (3^ fx, (2.13) 
that has the familiar time-varying coefficients 



a 



d f a' d' 



T{t) = - -2c + 4:d, a{t) =ab-cd + d^ + - (- . (2.14) 
a 2 \ a d 
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The time- dependent coefficients ao, 7o, ^o; ^o? satisfy the Riccati-type system fl2.6p - fl2.lip 
with Co = and are given as follows |13], |179] : 

«o(t) = v^#l-#„ (2.15) 

/3o W = A (t) = exp £ (c (s) - 2rf (s)) ds^ , (2.16) 

7o(t) = ^^ + ^ ° (2.17) 

2/ii (0) /io (t) 2a (0) ^ ^ 



and 



/io (t) 



/ is) - ^9 is)) /io (.) + #-/i'o (3) 



a (s) J 2a (s) 



ds 

W)' 



(2.18) 



2a (t) A (t) ^ , , ais) (J is) \(s) , , ^ ^ , , 

^o(t) = YTrf^'^o (t) + 8 / \\ ^ /io (^) <5o (.)) c^g (2.19) 

/^o l^J (/io (s)) 

io /^o(5) V a(s) y 



a (t) f^o (t) r2 u\ A f a(s)cT(s) 



«:o(i) = ^^^^TTTT^^o (i) - 4 / 77:,\7 (/io(^)^o(3))' (2.20) 
/^o io (/io (s)) 

""^"^ if^ois)So{s))(fis)-^gis)) ds 



/i^(s) V a(s)' 

(5o (0) = — eo (0) = g (0) / (2a (0)) and kq (0) = 0) provided that /ig and /i^ are the standard (real- 
valued) solutions of equation f l2.13p when cq = corresponding to the initial conditions /ig (0) = 0, 
/ig (0) = 2a (0) 7^ and (0) 7^ 0, fi[ (0) = 0. (Proofs of these facts are outlined in Refs. jl3] and 
|47] . Here, the integrals are treated in the most possible general way which may include stochastic 
calculus.) 

The systems f l2.6p - fl2.lip can be solved by the following variants of a nonlinear superposition 
principle |116j and [174J. 

Lemma 2. The solution of the Riccati-type system li2.6\) - l[KTl\) when cq = is given by 

/i(t) = 2/i(0)/io(t)(a(0) + 7o(i)), (2.21) 

/5(t)---4^^-^MA(t), (2.23) 
2(a(0)+7o(t)) ^{t) ^ ' 

'j(t)='j (0) - , , , ^ ' , (2.24) 
4(a(0) + 7o(t)) ^ ^ 



and 



Am A m /3o (t) (^(o) + ^o(t)) 

b [t) = do [t) --— , 2.25 

2 [a (0) +7o (t)) 
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'^'^ = 2(a(0)+7o(t)) ' ^'-''^ 

.(t) = -(0) + ^o(t)-^f^|^^ (2.27) 

4 (a (0) +7o (t)) 

m terms of the fundamental solution /I2.15\) - f2.2(j^) subject to the arbitrary initial data /i (0) , a (0) , 
/3(0)7^0, 7(0),5(0),£(0),«:(0). 

Lemma 3. The solution of the Ermakov-type system li2.6\) - l[2Jl\) when cq = 1 (7^ 0) is given by 

(2.28) 



/i = (0) /iov//3' (0) + 4(a(0) + 7o)', 

« = «o - Pi-, ^^^1±^ (2.29) 

'^°/3^(0)+4(a(0)+7o)' ^ 



/3 ^ ^ (t) , (2.30) 



(0) + 4 (a (0) + 7o)' ^ 



^ = ^(0)--arctan^^-^^^|^, a (0) > (2.31) 



^ ^ ^ _ . e (0) (0) + 2 (g (0) + 70) jS (0) + gp) ^ 
° /J^O) + 4 (a (0) + 7o)' 



25 (0) (« (0) + 7o) - /3 (0) (5 (0) + £o) 



(2.33) 



(0) + 4(a(0) + 7o)' 

^ = ^{0) + no-e (0) (0) ^ + "° (2.34) 

/3 (0) + 4 (a (0) + 7o) 

..^^ , e^O) (0) - (^ (0) + gp)^ 
+ (a (0) + 7nj -. ^ — 

/3nO) + 4(a(0)+7o)' 

m terms of the fundamental solution ^2.15\) - f2.20\) subject to the arbitrary initial data /i (0) , a (0) , 
/3(0) ^0,7(0), 5(0), £(0),«:(0). 

Here we would like to present a new compact form of these solutions. In order to do that, let us 
introduce the following complex-valued function: 



2a (0) + ^0 it) + ^ + (0) /io (t) (2.35) 

a(0)/ /ii(0) 

(a complex parametrization of Green's function and linear invariants of generalized harmonic oscil- 
lators are also discussed in Refs. [52] and [75]). Then 

z = CiE (t) + C2E* {t) , (2.36) 

where the complex-valued solutions are given by 

E{t) = ^^+ ^^^o (t) , E*{t) = l^ - ^/io (t) (2.37) 
/xi (0) /xi (0) 



6 ALEX MAHALOV AND SERGEI K. SUSLOV 

and the corresponding complex-valued parameters are defined as follows 

1 + f rf(0)\ 1-/3^(0) / rf(0)\ 

with 

Ci + C2 = l, |ci|'-|c2|' = ci-c; = /32(0). (2.39) 

In addition, 

2(0) = Ci + C2 = l, 2' (0) = 2za (0) (ci - C2) . (2.40) 
The inverses are given by 

|Cl| - |C2| |Ci| - \C2\ 

and 

/^i ^ (ci - C2) ^ + (ci - C2) 
2i{ci-ciy /ii(0) 2(ci-c^) ^ ■ ^ 

in terms of our complex vector (12.351) . 

One can readily verify that 

{cl-cl)z + {ci-C2)z* d(0) 
7o = KTT.. Z^^ + 



2i{z-z*) 2a (0) 

and equations (I2.18l) - (l2.20p can be rewritten in terms of vector z in view of (I2.42p . 

Finally, we introduce the second complex vector: 

C = £ (0) /3 (0) + « (5 (0) + £0) = C3 + ^£o, C3 = £ (0) /3 (0) + 2(5 (0) (2.44) 
and indicate the inverse relations between the essential, real and complex, parameters: 

« (0) = ^ - (0) = c, - c; = |c,r - \c.? (2.45) 



and 



5(0) = ^, .(0) = ±^^liS==. (2.46) 

2\/ |ci|^ - |C2|^ 



Then solutions of the initial value problem for the Riccati and Ermakov-type systems can be found 
in the following complex forms. 

Lemma 4. In terms of our time- dependent complex vectors ( and z, the solution of Riccati-type 
system is given by: // = yU (0) Rez = fi (0) {z + z*) /2 and 



Imz z- z* V |ci| - |c2| 

a = ao ^- — = ao 9 , (3 = ±2X^ , (2.47) 

2\zfRez 2i\zfz + z*' ^ z + z* ' ^ ' 

7 = 7(0)-:^^^, 5 = 5o + A^ = 5o-a^^, (2.48) 
2i z + z* tiez z + z* 



e = e(0)±i^/N'-N^^. '^^MOj + .o+gl^^. (2.49) 
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Lemma 5. The Ermakov-type system can be solved as follows: yU = /x (0) \z\ and 



a = ao + X'^^ P = ±X ^ "\ , ^ = ^(0)-larg;^, (2.50) 

2i\z\ z-z* \z\ 2 



S = So + X^^-^, e = ± + , (2.51) 

2^ \z\ o Ul /l^ |2 I |2 

II 2\Z\ \ Ci - C2 



[Cz+C z*) 



\z - z^ 



= (0) + /to + 7^ ,2 • (2.52) 

82 (ci - cl) \z\ 

(The proofs are rather straightforward and left to the reader.) 

As a consequence, for the Ermakov-type system one gets 

z + z* 





2i {a — ao) 






') + T 


z — z* 


.6 -So 






/3 (0)1^1' 


•■*{ 


P ) 



z — z" 



(2.53) 



i{a-ao) + '— = f—,, t{a-ao)-^ = P^^^ (2.54) 
..,-,* 2 z — z* 

,^0)Um±^Y (2.56) 



and 



-3(0) 

« = MO) + «„ + ^.-^!^t^.(0). (2.57) 

These "quasi-invariants" can be useful, for example, when making comparison of calculations done 
by different approximation methods. 

In this paper, we extend "The Simple Three Lemmas Approach", which is summarized and 
modified above for the linear problem, by the complexification of all time-dependent coefficients and 
by a classification of the simplest nonautonomous nonlinear terms. For a generic nonautonomous 
(derivative) nonlinear Schrodinger equation of the form 

i^t = Htp + R{tp) , (2.58) 

where H is an arbitrary variable quadratic Hamiltonian and R [tp) = P [tp, ip*, ip^, tpl) is a polyno- 
mial in four variables, we assume the natural gauge invariance condition: 

P (^e^^ re-'', i^e^'),, (^^e"^^) J = Ce^'P (^, r, C) • (2-59) 
The lowest terms that satisfy this condition are given by 

P (^, r, i^., O = hoi^ + ihx + h) + V'. + ^h,^p'^p: + h, l^l' ^ (2.60) 

where hk = {x, t) , = 0, 5 are some real or complex-valued functions. Our result is as follows. 

Theorem 1. The substitution ^2.4^ transforms the non- autonomous and inhomogeneous equation 
^2. 53^) with the lowest gauge invariant nonlinearities Ii2.60\) . namely, 

iip^ + a (t) ip^^ — h (t) x^ip + ic (t) xip^ + id{t)ip + f (t) xip — ig (t) ip^ (2-61) 
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into the autonomous form 

- iXr + X^^ - Co^^X = doX + idi^ + d2)\xfx + id3\x\'^X( + id4x'^{x^y + d::,\xfx 

(Co = 0, 1) 



(2.62) 



provided that 



hi = al3'^ 
h2 = |yu| 



die + d2 + -ds - -^"4 



2ImS 



h, = d,a^, l/ile^i-^ 
hf) = doa/S'^. 



(2.63) 
(2.64) 

(2.65) 
(2.66) 



Here do, di, d2, dz, d^, d^ are constants and S = ax^ + 5x + 



Proof. For the case of complex-valued coefficients, the transformation of the linear part is similar 
to [116] and |180j (with the use of contour integration if needed). Changing the variables in the 
nonlinear part: 

/i^/'e-'^P (2.67) 



h.X + ^^3|| e-^^-^ \X? X, + ./^4^e-^^- V {X,r + ^e'^^-^ 



x\ X 



+ 



h2 - hsS + hi5* + {hi - 2a/i3 + 2a*hi^ 



m 

-2ImS Ixl 



-X 



and substituting into fl2.6ip , one completes the proof with the aid of our conditions fl2.63p - p.66p . □ 

When cq = do = di = d^ = di = d^ = d2 = ho ^ and do = = ^2 = c^s = c?4 = 0, d^ = ho ^ 
0, we reproduce the results of Lemma 1 for p = 2 and p = 4, respectively. The derivative nonlinear 
Schrodinger equation, which was ffist derived for the propagation of circular polarized nonlinear 
Alfven waves in plasma physics [132] , |133] , and |192j , appears when co = do = di = d2 = d^ = and 
d^ = 2d/i (see [3], [2S], [21], [M], [SS], |14U] and the references therein for methods of solution 
of this equation). The amplitude equations of cubic-quintic type have been derived via asymptotic 
analysis of the governing Navier-Stokes equations of fluid mechanics in the limit of long spatial and 
slow temporal oscillations near the onset of instability [139j . Among other special cases are the 
Chen-Lee-Lui derivative nonlinear Schrodinger equation [26] and the Gerdjikov-Ivanov equation 
[65], [202] . 

Generic autonomous equations of the type 02.621) and some of their extensions such as quintic 
complex Ginzburg-Landau equation are discussed in [19], [22], [31], [32], [31], [60], [61], [62], [63] . 
[61], [83], [89], [99], [I03], [m], [Hi], [128], [I35], [136], [HI], [I62], [I63], [190], and [I93] (see also 
the references therein). 

If conditions fl2.63l) - fl2.66l) are not satisfied in certain application, one yet can find as functions 
of time and spatial variables thus moving the time dependence into the nonlinear part only. Then 
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one may use perturbation methods and/or some parameter control when possible. The Feshbach 
resonance in Bose-Einstein condensation provides a classical example of such nonlinearity control 
[50] . [88] . |152] . |153] . |154] . |156] . and |157] . A justification of the so-called local density approx- 
imation (see [19], |104] and the references therein) can be obtained from Lemma 1, cq = 0, under 
the following adiabatic condition: 

when a classical motion of the corresponding quadratic system is already taken into account (in 
general, the derivative should be taken with respect to a small parameter of the system under 
consideration). 

Moreover, in some important autonomous special cases, our Theorem 1 allows to determine the 
maximum symmetry groups of the particular equations; see, for example, |141j , |142j , |124j and the 
references therein. An extension to random- valued coefficients / and g is also possible, cf. [58] . 
Certain variations of initial data can be analyzed too. 



3. Integrability of Generalized Autonomous Nonlinear Schrodinger Equations 

In the case Cq = do = di = 0, a. detailed Painleve analysis of equation fl2.62p is performed by 
Clarkson and Cosgrove [53] (see also [22] for the extension to the case of complex parameters) . They 
have shown that this equation possesses the Painleve property for partial differential equations only 
if (is = ^4(2^4 — (i3)/4. When this relation holds, this is equivalent under a gauge transformation 
[TT3] . [89]: 

X = (j)exp (^-iu d^^ , (3.1) 

to a hybrid of the nonlinear Schrodinger equation and the derivative nonlinear Schrodinger equation, 

- z<p^ + + A 101' <j) + tfi (101' 0)^ = (3.2) 

(A, /i and u are constants), where one can assume that A = without loss of generality [193] . The 
corresponding derivative nonlinear Schrodinger equation is known to be completely integrable [3], 
[S], [2S], [M], [HS], [S5], [112], [HD], and [HD] (see also the references therein). Then solutions of the 
original equation are constructed by the gauge transformation (13.1 p in principle (see, for example, 
Refs. [33], [32], and [89] for more details). Explicit solitary wave solutions can be found in [59] , 
[89] . [97] . [99] . [117] . [128] . and [163] (see also the references therein). Our Theorem 1 allows to 
extend these results to a larger class of nonautonomous and inhomogeneous nonlinear Schrodinger 
equations. In the next sections, we will apply these general results to a nonautonomous solitons 
and other coherent structures in inhomogeneous medias including atmospheric plasmas. 



4. Cubic Nonlinear Schrodinger Equation 



The cubic nonlinear Schrodinger equation has a plane wave solution with an amplitude dependent 
phase. Under a simple condition (Lighthill's criterion), the modulation instability leads to self- 
focusing, namely, the parts of the wave front with the larger amplitude will propagate slower and 
this effect will enhance itself by the focusing from adjacent parts (Benjamin-Feir instability [18] . 
[166] . [177] , [182] , and [194] ). If the pressure of waves can balance the attraction force, the new 
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stationary solutions exist, which cannot be obtained by small nonlinear perturbation of the linear 
problem. These solitary waves, preserve their shapes during the time evolution. If the nonlinearity 
cannot balance the dispersion, a finite time collapse or dispersion that leads to decay occur. A 
detailed study of the exact solutions by different methods and many applications of the nonlinear 
Schr5dinger equations can be found in [2], [5], [30], [50], [63], [76], [#], [HB], [87], [I02], ^U\, [TT5] . 
[138] . [Us], [in], [M], and [180] (see also the references therein). 

Equation (12. 5 p is integrable when cq = and p = 2 |204j . |205j . and |206j . It has two standard 
forms, namely, focusing and defocusing. In this section, we concentrate on nonspreading and self- 
similar solitary wave solutions. 

4.1. A Similarity Reduction to the Second Painleve. In the traditional notations, the corre- 
sponding defocusing and focusing nonlinear Schrodinger equations, 

tij, + ij,, = ±2\ijfij (4.1) 

by the following substitution 

^ {x,t) = e^K-2^*V3)* gi/3p (^1/3 _ g,2^^ ^ ^ (4.2) 

(a is the acceleration) can be reduced to the (modified) second Painleve equations: 

F" = zF± 2F^ z = g^'^ [x - gt^) , (4.3) 

whose (bounded) solutions are the nonlinear Airy functions with known asymptotics as z ^ ±oo 
[33], [T71] . 

Combining with the familiar Galilei transformation |165j . 

^ {x, t) = e*(^-<2)^/2 ^ _ ^) (4_4) 

{v is the velocity), one obtains a more general solution of this type 

X g^/^F {g^^^ [x-vt- gt^) ) (4.5) 

in terms of the second Painleve Airy functions. (A similarity reduction of the nonlinear Schrodinger 
equation to the second Painleve equation was also discussed in Refs. [63|, [M], [171j . and [181]. The 
hnear case is investigated in Ref. [T7| where self-acceletating Airy beams were introduced.) With the 
help of Lemma 1 we extend these results to the most general nonautonomous and inhomogeneous 
integrable system of this kind. 

4.2. Self-Similar Solutions. The substitution 



vl/(X,T) = x( a/J-X,^T) (4.6) 



transforms equation (12. 5p into the focusing and defocusing forms in the new coordinates X and T : 

2^T + ^xx±2|^|^^ = 0. (4.7) 

By Lemma 1, these transformations give a similarity reduction of the nonautonomous nonlinear 
Schrodinger equation (12. 2p to the second Painleve equations (14.30 . As the result. 



^ (X, t) = _i_e<"W^^+^W^+''W) i ^t)x + e (t)) , ^7 (t) ] , (4.^ 
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where \l/ {X,T) is the solution, say in terms of the second Painleve transcendent in (14.51) . with a 
trivial change of the notation. 

As is known, the nonautonomous Schrodinger equation (12.21) under the integrability condition 
(12. 3p has also the following solution: 

ip {x,t) = —— exp (^i (^ax'^ + /3xy + ■y (^y'^ — go) + 6x + ey + nj) (4.9) 
V A'" 

xG {I3x + 2-fy + e) , 

where the elliptic function G satisfies equation 

dG \ ^ 1 

] = Co + 9oG^ + 2^oG"^ (Co is a constant of integration) . (4-10) 

and (j), y, Qq and Hq are real parameters (see also Ref. [174] for a direct derivation of this solution). 
Examples include bright and dark solitons, and Jacobi elliptic transcendental solutions for nonlinear 
wave profiles [2], |112] . |143] . |165] . |174] . In the original case (14. ip . setting Cq = y = 0, gives the 
stationary breather, which is located about a; = and oscillates at a frequency equal to go |164] . 
[T65] . 

4.3. Asymptotics and Connection Problems. In the defocusing case, the nonlinear Schrodinger 
equation (14. ip has the bounded solution, 

^P (X, t) = e*(--V2)./2+.,(.-.t-29iV3)t^l/3^^^ (^1/3 gt2^^^^ ^ (4.11) 

in terms of the nonlinear Airy function Ak^, when — 1 < /cq < 1 and ko ^ 0. The corresponding 
asymptotics are investigated in [5], [12], [33], [35], [37], [lO], [IH], [19], [130], [m], [158], [167] . 
and [183] (see also the references therein for study of this nonlinear Airy function). The end results 
are 

k()Ai (z) , z — 7- +00 
^'^oi^)={ r\z\-'/'sm(s(z)-eo)+o(\zr'/'), z ^ -oo. ^^"^^^ 



(the parameter represents amplitude of this nonlinear wave). Here, the following relations hold 

' ^ 2 ,^|3/2 _ 3^2 1^ ^ ^2 ^ 1^ (1 _ ^2) (413) 



and 



S [Z] = 

00 = ^r' In 2 + arg r A - ^r^") + ^ (1 - 2sign (fco)) • (4.14) 



If I/eqI = 1, then 



F (z) ~ sign (ko) v^^72, ^ ^ -oo. (4.15) 
If |A;o| > 1, then F (z) has a pole at a finite point z = Cq, dependent on ko, and 

F{z) resign {ko){z- Co)-', z^4. (4.16) 

In the focusing case, the reduction (14. 2 p results in the modified second Painleve equation, 

F" = zF- 2F^ (4.17) 
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and any nontrivial real solution satisfies |33j : 

F{z) =r\z\~^^'^ sin {s{z) -60) + 0(\z\~^^Un\z\y z ^ -00, (4.18) 

where 

s(z) = '^\zf^ + ^rHn\z\ (4.19) 
with r 7^ and arbitrary real constants. 
The second asymptotic is as follows. If 

^0 + ^r2ln2 - ivr - argr Qr^^ = 7™, n = 0, ±1, ±2, ... , (4.20) 

we have 

F (z) ~ /coAi (2) , z-^+00, (4.21) 



where ko is a nonzero real constant. The connection formulas are 

r" = TT 



2 - In (1 + kl) , sign (ko) = i-ir . (4.22) 



For the generic case, when the condition fl4.20p is not satisfied, one gets 

F (z) = a^/^ + a(3\2z\''^^^ cos {s{z) + 9) + 0{z-^), z ^ +00, (4.23) 

where a, (3 > and 9 are real constants, and 

s{z) = ^i2zf'-^pHnz. (4.24) 
The connection formulas for a, (3 and 9 are given by 

1 1 + If P 

a = -sign(ImO, /3' = -ln— (4.25) 

IT 2|lm4| 

= -l^- ^/3'ln2 + arg (l + e) + argT {z(3^) , 

where 

^ = (exp (vrr^) - l) exp Qr^ In 2 - ^tt + - arg T ^^^r^^ ^ ^ . (4.26) 

(See [33] and the references therein for more details.) 

In sections 7.1-7.2, we discuss applications of the second Painleve transcendents to nonlinear 
waves in fiber optics, plasmas and ocean. 

5. Derivative Nonlinear Schrodinger Equations 

The derivative nonlinear Schrodinger equation describes the propagation of circular polarized 
nonlinear Alfven waves in plasma physics |132] . |133] . and |192] . This equation is completely 
integrable [3], 0, [2S], [Si], ISH], [IIS], [I2D], and [HO] (see also the references therein). Its 
different forms are related through the gauge transformations |113] . |193] . Explicit solutions can be 
found, for example, in Refs. [59], [95], [96], [97], [99], [IH], and [T78] . 
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6. QuiNTic Complex Ginzburg-Landau Equation 



A number of nonintegrable nonlinear dissipative PDEs are known to display a wide variety of 
complex behavior, where the global time evolution is governed by the dynamics of spatially localized 
structures. This was shown in particular for the family of exact solutions of the one- dimensional 
supercritical complex Ginzburg-Landau equation [7], [H], [15], [21]. 

Exact solitary wave solutions of the one- dimensional quintic complex Ginzburg-Landau equation 
are obtained in Refs. [128] and |163] . In their notation and terminology, 

8 A d'^A 

— = eA + {h+ tci) Q^-ih- ^c^) \Af A - {h - tc,) \Af A, (6.1) 

where e, bi, Ci, 63, C3, 63, C3 are real constants and the field A{x,t) is complex- valued (see also 
[29j, [52], |172] . |173j ). These solutions are expressed in terms of hyperbolic functions, and include 
coherent structures with a strong spatial localization such as pulses and fronts, as well as, sources and 
sinks. Equation (16. ip is a one-dimensional model of the large-scale behavior of many nonequilibrium 
pattern-forming systems (see, for example, [7], [15], [21], [23], |128] . |162j . |163j and the references 
therein) . 

A systematic method for obtaining analytic solitary wave solutions of nonintegrable PDEs has 
been introduced by Conte and Musette ^40j, [137J and further developed by Hone [70], [80j and 
Vernov [W], [ISB] (see also [22] for another approach and [H], [12], [03], [100] for exact 

solutions). The unique elliptic traveling wave solution of (16. ip is found in Ref. [188j . Our approach 
allows to extend some of these results to nonautonomous systems under consideration. 



7. Some Applications 



For applications of Schrodinger equations in a variety of nonlinear problems, see Refs. [25], |28] . 

m, iza, m, nzi, m, m, m ra, ra, m\, 

[TB5] . [200], [102], [M], [^, and 

7.1. Accelerating Airy- Type Packets in Optics. In a Kerr medium with cubic nonlinearity, 
for which the dependence of the index of refraction on intensity is given by 

n = Uq (u) + ix (uj) + n2 \Ef , (7.1) 

the propagation of optical pulses in single mode dispersive fibers is described by the nonlinear 
Schrodinger equation of the form [25], [76], [77], |177] : 

^ + + .0^ J + 2-0^ + IV^I ^ = 0, (7.2) 

which is derived under assumptions that the complex envelope amplitude function ilj{x,t) varies 
slowly compared to the carrier and that the nonlinear and dispersion terms are weak. Here, 

and V is a geometric factor which depends on the radial variation of the guided electric field (we use 
the notation of Ref. [7^). In equation (17.21) . the second term describes the envelope propagation 
with group velocity uJq (in the absence of the rest of the terms), the third term describes the effect of 
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absorption, the fourth term describes the effect of dispersion, and the fifth term describes the effect 
of nonhnearity (further details can be found in [TB]). This equation is usually called a nonlinear 
parabolic equation in plasma physics |71], [HI], and |85] . 

In most of applications, equation (17. 2p has constant coefficients and, therefore, is not integrable 
if z/q 7^ 0. A similar special case of our generic equation (I2.2p . when a = uJq/2, b = c = f = 0, 
g = —oJq, and h = —vuon2/no can be transformed into the standard forms (14. ip by Lemma 1 with 
Co = 0. In our notations, the standard solutions are given by /Zq = 2ate'^'^^ and /x^ = (1 — 2dt) e^*^*. 
By fl27[5|) - (l2:20|) . A = e^^* and 

«o (t) = -2/3o (t) = 7o W = ^> 5o = -^o = £, ^o{t) = ^. (7.4) 
Then /x = /i (0) (1 + 4a (0) at) e^^^ and 

^(°) ^^^(0)- ^'^W"' (75) 

l + 4a(0)at' l + 4a(0)at' l + 4a(0)at' 



1 + 4a (0) at ' ' ' ' ' 1 + 4a (0) at ' ' ' ' ' 1 + 4a (0) at 

by Lemma 2, thus providing the general transformation (14.80 into the standard forms (14. 7p . Finally, 
the integrability condition (12. 3 p takes the form 

h = hof3^ (0) ^ (0) = a(l + 2t(d-2a(0)a)t) + (t^) , t ^ 0, (7.6) 

1 + 4a (0) at 

which indicates the most general integrable system of this kind and may be achieved in fibers, 
say, through the geometrical factor v in order to support nonspreading signals. Moreover, our 
equations (17. 5p describe explicit evolution of variations of initial data. The choice of initial condition 
d = 2a (0) a allows to include effects of absorption into the zero approximation. Our modification 
of initial condition allows to compensate some intensity losses. 

The bright and dark soliton solutions are given in Refs. [76] and [77] (see also [21], [136] . [204] - 
and the references therein). But here, we mainly concentrate on solutions related to the second 



Painleve transcendents. 

Explicit solutions (14. 2 p and (14. 5 p in terms of the nonlinear Airy functions reveal several remark- 
able features. In quantum mechanical terms, the "probability density" \ip {x,t)f not only remains 
unchanged in form but also continually accelerates in empty space. Therefore, these solutions ex- 
tend to nonlinear cases the nonspreading Airy packet introduced by Berry and Balazs [17] for the 
unidimensional linear Schrodinger equation without potential (see also [72] and [186] ). There is 
no violation of Ehrenfest's theorem because the Airy function is not square integrable and cannot 
represent the probability density for a single particle. What accelerates in the Airy packet is not 
any individual particle but the caustic (i.e., the envelope, or focus) of the family of orbits (see [IT] 
for more details on this geometrical approach). These nonspreading and freely accelerating wave 
packets have recently been demonstrated in both one- and two-dimensional configurations in optics 
|169] . [170] and we would like to elaborate on this topic in detail. 

It is worth noting that an expansion transformation from the Schrodinger group, e.g. formula 
(2.8) of [124] with m = — 1/ti, gives another accelerating solution of the free particle equation, 
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iip^ + ip^^ = 0, in terms of Airy function as follows: 



/ \ti\ (. ( 2g tie \ tit ix" , 



xg'/'Ai{g'/^x-g 



tit'' \ ti 



ti-tj ti-t^ 

which is convenient for t < ti. The degenerate case, when ti = 0, can be analyzed with the help of 
transformation (2.9) of |124] : 



From now on, we choose ti > for the sake of simplicity. (The most general solution of this kind 
can be obtained by transformation (2.5) of |124] .) 

An arbitrary point Ai(xo = constant) on the Airy function graph accelerates in empty space 
according to the law: 

^ = ^ - ^9tl + (^gti + ^) (ti - t) , (7.9) 
dx gt\ / Xq 

(^7^"r' + ^ 

rf^x _ 2gt\ 
^ ~ {ti - tf 

in the laboratory frame of reference. This "kinematics" can be obtained as a unique solution of the 
following singular boundary value problem on (— cxD,ti): 

(^1 - rf^ - ~ ^) - = (7.10) 
(an analog of Newton's second law of motion) , 

lim P^,=Ci = gti + ^ (7.11) 

t^-oo ti — t g ' ti 

(dx 
— ~ — Ci, an asymptotic value of constant velocity at — oo 
dt 



\im {ti - t) X {t) = C2 = gtl ^ (7.12) 



(a; ~ C2/ {ti — t) , a given residue of the simple pole at ti) . 
The corresponding family of orbits in phase space is given by the equation: 

^ = 2Ci + P + ^, (7.13) 

where P = dx/dt + Ci and Q = x + 2gt\. By the quadratic formula. 
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This curve degenerates to a parabola when Ci = 0. 

In view of fl7.9p . any two points on the graph, say Ai(a;o) and Ai(?/o) with xq < Hq taken at the 
same initial moment of time to? are contracting to each other: 



x{i)-y (t) = - yo) 



U - 1 







as t — )■ ti 



(7.15) 



and their trajectories have two different slopes as t — ?■ ±oo, when dispersion that leads to decay 
occurs. In other words, these analogs of classical trajectories have slanted asymptotes and linearly 
diverge from each other as t — )■ ±oo : 



lim 



9^0 + ^] {ti -t) + 2gti 



0, 



Mm {{ti - t) X (t) - gtf) = (7.16) 



but both of them approach to a common vertical asymptote at t — )■ ti. With the aid of Airy function 
asymptotic [146] . for any fixed spatial point x one gets 



\^P{x,t)\r.gy^^[(t,-tfz\ 



-1/4 



sm -z 
3 



3/2 



as t — 7- > 0, 



where 



9 



1/3 



gUf 



ti 



> 0, 



lim (ti 



tfz 



g''H\. 



(7.17) 



(7.18) 



ti-t J h-t 

Thus the leading term remains bounded and does not depend on the spatial coordinate x when 
t — 7- . As a result, a possible "blow up" of this solution cannot occur at any finite point in space 
as expected in the linear case. 

On the contrary, the cubic nonlinear Schrodinger equation is no longer preserved under the 
expansion transformation (2.8) of [124J. But the same symmetry, say by our Theorem 1, holds for 
the quintic nonlinear Schrodinger equation, which is thus invariant under the action of this group 
of transformations. This is where the blow up solutions do exist and some of them are analyzed in 
section 7.5 (see, for example, equation ([77 



Our solution (17. 7p can, in principle, be generated from the Airy packet of Berry and Balazs [T7] 
in the following fashion. Let us consider a generalized harmonic oscillator with the variable Hamil- 
tonian (12.21) . which evolves from the initial free particle coefficients and terminates its evolution to 
the original case of a free particle. Solution of the corresponding Riccati-type system is provided 
by Lemma 1. By the end of the cycle, the original Airy packet may possess the initial data of (17. 7p 
and this solution will be engaged by the continuity of the wave function. In addition, we provide an 
example of instability for the linear Airy beam, namely, the corresponding variation of the initial 
configuration: 



|ti 



ti — to 



i/^Ai 



X exp \ tg \ X 



2^7 



,1/3 



titl 



X 



ti — to J ti — to 



(7.19) 



tfto 



IX 



3 ti-toj {ti-toY 4(ti-to). 

results in a "collapse" of the packet in a finite time to < t < ti. In this process, any point from 
the original beam attains an infinite velocity and acceleration as t — )■ in the laboratory frame of 
reference due to deterministic evolution of our solution (one may call this exotic dynamic an "Airy 
gun" ; a Mathematica animation is available on the article website) . 
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The study of propagating nonlinear Airy-Painleve optical pulses in dispersive fibers was initi- 
ated in [nn] (see also Ref. |171] for an earlier application of the second Painleve transcendent in 
hydrodynamics) and has been continued in recent publications [I3], |90], |126] . and |161] . Our 
solution (14. 8p incorporates the simplest nonhomogeneous effects of a nonlinear time-varying media 
in a unified form. The corresponding asymptotics, connection problems, and bibliography are given 
in section 4.3 for the reader's convenience. This summary may facilitate further use of these results. 
An important example of motion of the nonlinear Airy packet in a time-varying spatially uniform 
force is left to the reader (see \V7\ for discussion of the classical case). 

It is worth addressing a few new features of the nonlinear Airy beams under consideration. 
According to (14.121) . in the defocusing case, which is related to positive group velocity dispersion 
(dark pulse), a bounded on the entire real line solution corresponds to the real parameter \ko\ < 1, 
/cq 7^ (a finite range of the ratio of the nonlinearity to dispersion [BS])- In the focusing case 
(anomalous dispersion, bright pulse), the bounded solution can only exist under condition f l4.2Up . 
which may be thought of as a "quantization rule" in this nonlinear problem. The connection relation 
take the form 



dr. 



ivr - ^ In (1 + kl) + arg r ( ^ In (1 + k^) ) + 



4 



271 



7m, 



ko^O 



(7.20) 



and the asymptotics of bounded solutions of the modified second Painleve equation are given by 



F(z) 



ln(l + fco') 



vr 



-1/4 



koM (z) 



sin (s (z) — 9r 



z — )■ +00, 



O 



z\ ^''^InU 



z — )■ — oo 



(7.21) 



s(z 



s{z) 



(7.22) 



with 

2 I |3/2 3^ 

3 47r 

for n = 0, ±1, ±2, .... For the bright pulse, there is no restriction on the nonlinear wave amplitude 
and this case deserves an experimental study. 



\zf^-—\n{l + k'^) \n\z\ 



In our opinion, a joint consideration of (breaking/restoring) symmetry of the linear and (cu- 
bic/quintic) nonlinear Schrodinger equations provides a natural hierarchy of accelerating, non- 
spreading, and blowing up wave packets. 



7.2. Nonlinear Airy Waves in Plasma and Ocean. In plasma physics, evolution of weakly 
nonlinear, strongly dispersive, short wavelength electron Langmuir waves described by their modu- 
lating envelope is governed by the Schrodinger equation with cubic nonlinearity. The corresponding 
solitons are called Langmuir solitons ^13"5] . |143j in the case of a homogeneous media. In laser 
plasma experiments, the plasma is both inhomogeneous and nonlinear to the electromagnetic waves 
and the large-amplitude plasma waves (see, for example, [27] and the references therein for more 
details). We demonstrate our approach in the simplest inhomogeneous situation. 

Solitons moving with acceleration in linearly inhomogeneous plasma were studied in Refs. 
and [28] (see also [9J and |184] ). Single solitons and multisolitons are found accelerated whereas 
maintaining their shapes when moving around even upon emerging from collisions with other soli- 
tons. 
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Here we analyze another interesting "one soliton" scenario that is somewhat complementary to 
recently discovered Airy beams in optics. In the defocusing case, 



^l^ + ^ = m'x- (7.25) 



the following simple substitution 

^P{x,t)=e-^<''^^^'"''l')x{i.T), i = x + 2kt\ T = t, (7.24) 
which is originally due to Tappert [27] (see also |180j for a more general case), results in 

When g = 2k, the composition of this transformation with solution (14.111) predicts that, despite of 
the presence of the external constant force, the following wave packet: 

{X, t) = ei(--t/2)-/2-.5^'*V2^1/3^^^ _ (7 26) 

moves with the constant speed v (the standing wave solution, originally found in Ref. |171] . occurs 
if w = 0). Here, the soliton profiles are defined as the Airy-type solutions of the second Painleve 
equations and the corresponding asymptotic properties of these nonlinear Airy functions are dis- 
cussed in section 4 for the both, defocusing and focusing, cases (see also [33] for more details). It 
is worth noting once again that in the focusing case the "quantization rule" f l4.20p should hold for 
the bounded solution. 

For an application of the second Painleve transcendents, in the case of the nonlinear Schrodinger 
equation with the linear potential (I7.23p . to giant waves, as observed on the Agulhas Current of the 
southwest Indian Ocean, see Ref. |171] . The operation of an electrostatic spherical probe in a slightly 
ionized, collision-dominated plasma can also be described by the second Painleve transcendents |3B] , 
|92] . and [03]. A similar behavior occurs in numerical simulations of the shock waves in plasma and 
corresponding experiments (see [SI], [SJ], [HI] and the references therein). 

7.3. Nonlinear Alfven waves. The propagation of a nonlinear Alfven wave with a small but 
not-vanishing wave number along the magnetic field in cold plasmas is governed by the deriva- 
tive nonlinear Schrodinger equation |132j . |133j . |192j . |191j . The Alfven solitons are discussed in 
|134j . |155j . |200j . The derivative nonlinear Schrodinger equation is also used for modeling of wave 
processes in nonlinear optics. Stokes waves in fluids of finite depth, etc. Our approach give an 
opportunity to incorporate the inhomogeneous effects of the media. 

Freak waves in plasmas and optical rogue waves are discussed in recent publications [S], |118] . 
[T60] . [T68] . and [1^ . 

7.4. Ion Acoustic Envelope Solitons in Explosive Ionospheric Experiments. The theory 
of nonlinear wave modulation in cohesionless plasma on the basis of the Vlasov description |189] 
with the nonlinear Landau dumping is developed in [81] and [82] (see also |196] ). The corresponding 
nonlinear Schrodinger equation that describes the amplitude of the plasma density fluctuations has 
a nonlocal- nonlinear term [82] . 

The stable ion acoustic envelop soliton propagating perpendicular to the magnetic field lines can 
exist in the ionospheric plasma. These solitons were identified by processing the data from the North 
Star active explosive ionospheric experiment [109] , |110] . The parameters of the soliton have been 
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first estimated under assumption that the coefficients in the Schrodinger equation are constant [109] . 
The additional ion dissipating terms have been used in order to obtain the stable sohton solutions 
in a plasma with steep density and temperature gradients |110j (see also the references therein). 

Our approach give an alternative opportunity to incorporate nonautonomous dissipative terms. 
Exact solutions of the cubic complex Ginzburg-Landau equation are to be taken as a nonpeturbed 
approximation. These solutions include coherent structures with a strong spatial localization such 
as pulses and fronts, as well as, sources and sinks [7], [39], [H], [15], |144] . |145] . |150] . and |163] . 
More details will be discussed elsewhere. 



7.5. Oscillating Coherent Structures. Let us consider the autonomous generalized nonlinear 
Schrodinger equation fl2.62p and address an intriguing question — can self-oscillating solutions 
exist? By our Theorem 1, the following special case Cq = 1, do = di = d2 = and d^ = d/i of this 
equation, namely, 

+ - x^ij = tds V'. + ^V:) + 4 1^1' ^, (7.27) 
(we replace r = —7 and use the original notation for convenience) is invariant under the action of 
the Schrodinger group, which was originally introduced by Niederer [142] for the linear harmonic 
oscillator, when d^ = d^ = and space-oscillating solutions exist |125] . |127] . Therefore the non- 
linear equation may possess some interesting oscillating solutions too! The direct action of the 
Schrodinger group is given by 



^ (x, t) = ^/H^ ,^iait)^^+sit).^m) ^(p[t)x + e it) , -7 it)) 



(7.28) 



where, according to Lemma 5, we define z (t) = cie^** -|- C2e and find everything in terms of this 
complex- valued function as follows 



"W = ^ '^'^''"":;:',f' ■ Ht)-± ^ O^W = iarg.. (7.29) 

Z yCl — €2) \z\ Z \z\ ^ 

c^z-c%z* ,^ c^z + c%z* iclz+dfz*){z-z*) 
Ht)= ' ., I , e{t) = ± ' ' K{t)= ,1 ,2 ■ 



2 \z\ y |Ci| - |C2 

The complex parameters: 

1 + (0) 1 - CO) 
c, = —^j^-ta{0), C2 = ^^ + ta{0), C3 = e (0) /3 (0) + z5 (0) (7.30) 

are defined in terms of 'essential' real initial data (here we choose 7 (0) = k (0) = for the sake of 
simplicity). The real form of transformation fl7.28p and visualization of the corresponding oscillating 
"missing" solutions for the linear harmonic oscillator can be found in Ref. [125] . In view of 

\zf = \c^f + \c^f + ci4e^'' + clc2e-^'\ cic; = ^—^^ - {0) - la (0) , (7.31) 

the shape-preserving configurations, that are somewhat similar to the coherent states of the linear 
harmonic oscillator, may occur through this transformation only when a (0) = and /3 (0) = ±1. 

Our formulas f l7.28p -( 17.30p provide a new time-periodic solution of equation f l7.27p from any given 
solution. Although explicit solutions of (I7.27P are not readily available in the literature (see, for 
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example, [38] and the references therein for ^3 = 0), Bose condensation and/or nonhnear effects in 
"non-Kerr materials", e.g. fiber optics beyond the cubic nonlinearity, provide important examples. 

7.5.1. Example 1. Mean-field theory has been remarkable successful at description both static and 
dynamic properties of Bose- Einstein condensates [50] , |154] . The macroscopic wave function obeys 
a 3Z) cubic nonlinear Schrodinger equation, which is usually called the Gross-Pitaevskii equation in 
this model. There are several reasons to consider higher order nonlinearity in the Gross-Pitaevskii 
equation [154j . The quintic case is of particular importance because near Feshbach resonance one 
may turn the scattering length to zero when the dominant interaction among atoms is due to three- 
body effects (see [19], |103j . |136] . |153] . |156j . |157j and the references therein; on Li, for example, 
the scattering length is reported as small as 0.01 Bohr radii |157j). Then the nonlinear term in the 
mean-field equation has the quintic form. Another interesting example is a ID Bose gas in the limit 
of impenetrable particles [68], |104] . |185j . 

The ID-quintic nonlinear Schrodinger equation without potential in dimensionless units, 

iAt + A^^±-^\A\'' A = Q, (7.32) 

has the following explicit solutions adapted from |128] (we use the notation and terminology from 
[128] and [163]; see also [60] and [63]). 



Pulses: 

1/2 



A (x, t) = e 



i(j> 



X 



cosh k{x — vt) 
exp^(^+(A;2-t;2) 



(7.33) 



(0, V and k are real parameters, the upper sign of the nonlinear term should be taken in f l7.32P ; see 
also |149] and the references therein). We have 



A{x,t)\^ dx = 7i (7.34) 



and the corresponding plane wave expansion 

1 

^271 



A (x, t) = -L / e'P'^B {p, t) dp, (7.35) 



can be found in terms of gamma functions: 



B [p, t) = -= exp i I pv ] t (7.36) 



2nVk V 4 



4 2k V 2) V4 2A; 2 



with the help of integral (1A.3|1 from Appendix (the case = allows to evaluate the L^-norm of 



solution fl7.33p ). The energy functional is given by 

^ = -V?-^H^=^>0 (7.37) 

and its positivity provides a sufficient condition for developing of a blow up, namely a singularity 
such that the wave amplitude tends to infinity in finite time [23], |177] . and [208]. 
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Sources and sinks: 

cosh (v^r (x — vt)) ± 1 



A{x,t) = e'^r 



1/2 

(7.38) 



cosh (v^r {x — vtj) =1= 2 
xexpi (^y - (t;2 + 3r2) 

(0, f and r are real parameters; see also [104] ). Equation f l7.32p has also a class of (double) periodic 
solutions in terms of elliptic functions [60], [63]. They will be discussed elsewhere. 

In addition, direct action of the Schrodinger group |124] . |141] on fl7.33p produces a six-parameter 
family of square integrable solutions: 



l + 4a(0)t V l + 4«(0)t 

(one can choose v = withot loss of generality; see also [6], [70], [73] and the references therein 
regarding the Schrodinger group). All of them blow up at the point Xq = —6 (0) /2a (0) in a finite 
time, when t — )■ to = ~l/4a (0) and a (0) ^ 0. (We use real- valued initial data for the corresponding 
Riccati-type system; see [124J for more details.) Blow up of these solutions can be naturally studied 
in phase space. The corresponding Wigner function [78], [197] is evaluatited with the help of integral 

(jAU). 

A finite time blow up of solutions of the unidimensional quintic nonlinear Schrodinger equation 
f l7.32p is a classical result discussed in many publications; see, for example, [20], [23], |149] . |159] . 
[207j . [208] and the references therein. This case is critical because any decrease of the power of 
nonlinearity results in the global existence of solutions [67], |177] . |195j . |208] (see also Refs. |105] 
and |104j for a renormalization approach). Here, for the family of solutions (I7.39p . we have shown 
that the blow up occurs due to a "hidden symmetry" of this nonlinear PDE. This property holds 
for all solutions which exponentially decay at infinity. 

7.5.2. Example 2. The quintic nonlinear Schrodinger equation in a parabolic confinement, 

ti;, + ^,,-x'i;±^\i;\'i; = 0, (7.40) 

describes a mean-field model of strongly interacting ID Bose gases for the practically important case 
of a harmonic trap [19], [53], [69], [104], [105] , [153] . [T59] and, in particular the so-called Tonks- 
Girardeau gas of impenetrable bosons [68j, [185] : sec [101] and [148] for experimental observations. 
(The time-independent version of the quintic nonlinear Schrodinger equation has been rigorously 
derived from the many-body problem [122] : see also [121] for a rigorous derivation of the Gross- 
Pitaevskii energy functional.) 

Our observation is as follows. By the gauge transformation (e. g. [B], [75], [124] and the references 
therein for the linear problem, the quintic nonlinearity is invariant under this transformation by 
our Theorem 1), equation (I7.40p has the following solution: 

-(^/2)xH^n2t / tan2t\ 

^ x, t) = A -, ^— , 7.41 

Vcos2t Vcos2t 2 J 
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where A {x,t) is any solution of (I7.32p . in particular, the pulses and sources (I7.33p and fl7.38p . 
Oscillating pulses: 

2vx + (fc^ — f ^ — x^) sin 2t 

X exp i 

^ 2cos2t 

(0, V and A; are real parameters, the upper sign should be taken in the nonlinear term). They are 
square integrable at all times: 

oo 

\^ix,t)f dx = 1 (7.43) 



1 

TT 

and 

- / \x'^(x,t)f dx = ^ cos^ 2t + sin^ 2t, (7.44) 

TT y_oo (4A;) 

- J_j,pj.,t)\^ = _,sm^2*+„^cos'^2« + ^^. (7.45) 

The corresponding elementary integrals are collected in Appendix for convenience; see (]A.4p . 
The expectation values and variances of the position x and momentum p = i~^d/dx operators 



are given by 

__ < 

— ' - (1) 



X = ^ = vsm2t, p= = vcos2t (7.46) 



and 



^2 ^2 

(Sx? = x^-x^ = ^cos2 2t, (Sp?=p'^-p^ = ^sin2 2t+ —, (7.47) 

^ ' (Ak)^ ' ^ ^ ^ (4^)2 2cos2 2t' ^ ^ 



x = t;t, p=^, (5x)^ = ^, {Spf = ^ (7.48) 



respectively. [It is worth noting that 

^ 2' ^ 4fc2' ^ 8 

with 

{6pf {6xf = ^ > ^ (7.49) 
for the original traveling wave solution fl7.33p ]. The energy functional is given by 

1 ^2 

E = H = p^ + x'^ - = ^ + w2>0 (7.50) 

A {Akf 

by the direct evaluation. 

A remarkable feature of the oscillating solution (I7.42p is that the corresponding probability density 
converges, say as a sequence, periodically in time, to the Dirac delta function at the turning points: 
\'(p (x, t) 1^ — )■ 7r5 (x =F f ) as t — )■ ±7r/4 etc., when an "absolute squeezing" and/or localization, namely 
5x = 0, occurs with 6p = oo. The fundamental Heisenberg uncertainty principle holds 

(M^(fa)^ = g(l + 3gl-^«)>g>l (7.51) 
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at all times. (It is worth noting that 7r^/8 ~ 1.2337. The minimum- uncertainty squeezed states 
for a linear harmonic oscillator, when the absolute minimum of the product can be achieved, are 
discussed in Ref. 

The corresponding Wigner function [78], |197] : 

1 

r {x + y/2,t)ij{x-y/2,t)e''Py dy, 



W {x,p,t) 



2ti 



can be evaluated in terms of hypergeometric function: 

, 1/2 + 1/2 

W {x,p, t) = sech CO 2F1 



— sinh d 



(7.52) 



(7.53) 



with the aid of integral flA.ip . Here 

1 



{p cos 2t + X sin 2t 



2k 



{x — V sin 2t) . 



2k ' ' ' cos2t 

A visualization of Wigner's function in phase space can be found at the article website. 



(7.54) 



Our mathematical example reveals a surprising result that a medium described by the quintic 
nonlinear Schrodinger equation (17.401) may allow, in principle, to measure the coordinate of a 
"particle" with any accuracy, below the so-called vacuum noise level and without violation of the 
Heisenberg uncertainty relation, which is a major obstacle, for example, in the direct detection of 
gravitational waves [T], [55] . 



Oscillating sources and sinks: 



^{x,t) 



2r 



31/2 cos 2t 
2vx 



cosh 



2r 



cos2t 

^2 _|_ ^2 _|_ gj^2t 



[X 



V sin 2t)) +2 



1/2 



(7.55) 



X exp t 

^ 2 cos 2t 

(0, V and r are real parameters, we have chosen the lower sign of the nonlinear term in f l7.40p ). 
Their detailed investigation will be given elsewhere. 

Then the action of Schrodinger group, say in our complex form f l7.28p -( l7.30p . on ( I7.42p and/or 
f l7.55p produces a six-parameter family of new oscillating solutions of equation f l7.40p . For example, 
the following extension of fl7.42p holds: 



ip{x,t) 




2k(3 (0) 



2a (0) sin2t + cos2t 



(7.56) 



,1/2 , I r,^ , X - (5 (0) + 17/3 (0)) sin 2t 



X exp I 



X exp I 



2a (0) sin2t + cos2t 
(2a (0) cos 2t - sin 2t) x"^ + 5 (0) {2x - 5 (0) sin 2t) 
2 (2a (0) sin 2t + cos2t) 
2v{x-S (0) sin 2t) + {k^ - v^) (5 (0) sin 2t 
2 (2a (0) sin 2t + cos2t) 



/3(0) 



+ ve (0) 
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which presents the most general solution of this kind. (We assume that 7 (0) = k (0) = for the 
sake of simplicity. Although the breather /pulsing solution, when a (0) = 5 (0) = e (0) = f = and 
(3 (0) = 1, was found in Ref. |159] . our discussion of the uncertainty relation and Wigner function 
seems to be missing in the available literature.) The blows up occur periodically in time at the 
points 

(5(0) +f/3(0) , , , , , 

xo = ± ] ' =, when cot 2t = -2a . 7.57 
v/4a2 (0) + !' 

The corresponding Wigner function is given by our formula fl7.53p with the following values of 
parameters: 

00 = I, . [{p - 2a (0) x) cos 2t + {2a (0) p + x) sin 2t - 6 (0) - v(3 (0)] , (7.58) 

V 2a(0)sin2t + cos2t ^ ' 

(One can put t> = in fl7.56p - fl7.58p without loss of generality because the general action of the 
Schrodinger group already includes the Galilei transformation |141j . |142j .) 



7.6. A Nonlinear Harmonic Oscillator and Painleve IV. The following special case of equa- 
tion ^rm -. 

i^t + V-xx - = 2x + 3 V, (7.59) 

by the substitution 

^(x,t) = e-^("+^/2)t^(x) (7.60) 

can be reduced to the ordinary differential equation for a nonlinear harmonic oscillator studied in 
Refs. [lO], dU, and [33]: 

u" = 3u^ + 2xu^ + -n-^^u. (7.61) 

There are exact 'bounded state' solutions for any integer n = 0,1,2,... that asymptotically ap- 
proach the wave functions of the linear harmonic oscillator as |x| — )■ 00 (see |33] for a summary 
of these results). The probability density \tp{x,t)f = u"^ (x) can be found in terms of a special 
fourth Painleve transcendent. Equation fl7.59p arises also as a symmetry reduction of the derivative 
nonlinear Schrodinger equation, which is solvable by inverse scattering techniques [3], (TU], and 
[9^ . (A symmetry reduction of the cubic-quintic nonlinear Schrodinger equation to fourth Painleve 
transcendents is discussed in p3].) There are a very few nonlinear ordinary differential equations 
which obey the Painleve property and have bounded solutions and this is one of those. It would be 
interesting to find a connection of these solutions investigated in Refs. [10]-[TT] in details with the 
theory of Bose-Einstein condensation and fiber optics, say for a stable signal propagation. 
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Appendix A. Integral Evaluations 



The following integral 

e'""' ds 



f 

J — c 



TT 



Vcosh s + cosh c cosh vrw 



1 1 



sinh 



sinh 



< 1, 



(A.l) 



which Mathematica fails to evaluate in a compact fornQ, can be derived as a special case of integral 
representation (2) on page 82 of Ref. [57]. The hypergeometric function is related to the Legendre 
associated functions, which are a special case of Jacobi functions, see [51], |106] . |107] : 



Pi/2-iuj (coshc) = 2F1 



sinh - 



(A.2) 



and Mehler conical functions. An important special case is given by 

T{l/4: + iuj/2)T{l/4:-iuj/2) 



00 gJOJS 



, a/ cosh s 
Useful elementary integrals, 

du f°° u^ du 



ds 



cosh-u 



vr, 



cosh u 4 



27r 



sinh u du Svr 



cosh u 



du 



cosh u 



IT 



can be verified by Mathematica. 



(A.3) 



(A.4) 
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